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ABSTRACT 

In this paper, the third of a series, we study the growth rate and luminosity of black holes (BHs) in motion with 
respect to their surrounding medium. We run a large set of two-dimensional (2D) axis-symmetric simulations 
to explore a large parameter space of initial conditions and formulate an analytical model for the accretion. 
Contrary to the case without radiation feedback, we find that the accretion rate increases with increasing BH 
velocity Vhh reaching a maximum value at v^h = 2c s .i n ~ 50 km/s, where c Si i n is the sound speed inside the 
"cometary-shaped" H II region around the BH, before decreasing as v^. The increase of the accretion rate with 
i?bh is produced by the formation of a D-type (density) ionization front (I-front) preceded by a standing bow- 
shock that reduces the downstream gas velocity to transonic values. Since the I-front is beyond the classical 
Bondi radius for the hot ionized gas, the accretion flow in the BH frame of reference is similar to the stationary 
case. Interestingly, there is a range of densities and velocities in which the dense shell downstream of the bow- 
shock is unstable; its central part is destroyed and reformed periodically, producing a periodic accretion rate 
with peak values about 10 times the mean. This effect can significantly increase the detectability of accreting 
intermediate mass BHs from the interstellar medium (ISM) in nearby galaxies. For «bh > 2c s in , the central 
part of the bow-shock is not able to regenerate, the I-front becomes R-type and the accretion rate approaches the 
classical Bondi-Hoyle-Lyttleton solution. We find that the maximum accretion rate for a moving BH is larger 
than that of a stationary BH of the same mass, accreting from the same medium if the medium temperature is 
< 10 4 K. This result could have an important impact on our understanding of the growth of seed BHs in 
the multi-phase medium of the first galaxies and for building and early X-ray background that may affect the 
formation of the first galaxies and the reionization process. 

Subject headings: accretion, accretion disks - black hole physics - dark ages, reionization, first stars - hydro- 
dynamics - methods: numerical - radiative transfer 



1. INTRODUCTION 

Gravitationally driven gas inflow onto moving point masses 
such as black holes (BHs) or neutron stars has been described 
analytically in the 40s by Bondi-Hoyle-Lyttleton ( |Hoyle & 
|Lyttletonll93"9l|Bondi _& Hoyle|1944[|Bondi|1952| l. The gen- 
eralized formula for the accretion onto a point mass moving 
with velocity Vhh is obtained from the Bondi formula by re- 



placing the gas sound speed, 



with an effective speed 



v e & = ( c s,oo + v bh) ■ The accretion rate M oc pv c s<r c tt is 
roughly estimated as the gas flux through the effective cross 
section er e gf = 7rr^ ff , where r e g = GMbh/ v eS i s trie distance 
at which the escape velocity of the gas equals v e g. Despite 
the similarities between the Bondi and Lyttleton formulae, in 
the second case the accretion onto the BH is not spherically 
symmetric: most of the gas streams past the BH and is grav- 
itationally focused on the axis of symmetry of the problem. 
The component of the gas velocity perpendicular to the direc- 
tion of the BH motion is converted into thermal energy and is 
mostly dissipated, thus a fraction of the gas becomes gravita- 
tionally bound to the BH and is accreted from the downstream 
direction. The generalized Bondi-Hoyle-Lyttleton formula is 
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where Mb is the Bondi accretion rate onto non-moving BHs 
Mb = Ke 3 / 2 pozG 2 M^cJ^z (assuming isothermal equation 
of state: 7 = 1). Here, Mbh is the BH mass, and poo is 
the density of the ambient gas. The term in the denominator 
of Equation ([Tji is the only term that accounts for the motion 
of the point mass, thus for supersonic BH motion the accre- 
tion rate decreases with increasing velocity as v^. Numeri- 
cal simulations of accretion onto moving BHs have confirmed 
the validity of the Lyttleton equatio n ( Shima et al.|1985 Ruf-| 
fert & Arriett||1994| |Ruffert|[P996] l, even in the presence of 



non-axisymmetric flows that necessarily arise in 3D simula- 
tions due to hydrodynamical instabilities (|Cowie|1977)|Mat-| 
|suda etal.|1987||Fryxell & Taam|1988||Taam & Fryxel l|i988| 
|Soker|1990||Koide et al.|1991||Livio et al.|1991[|Fogrizzo &| 
Ruff ertl 1 997| 1 1 9^>9j |Foglizzo et al.|2005f ~^ 

For the case of accretion onto stationary BHs, the radia- 
tion emitted near the gravitational radius of the BH typically 
reduces the rate of gas supply to the BH from large scales 
well below the value given by the Bondi formula. This is be- 
cause X-ray and UV heating increases the sound speed of the 



gas in th e proximity of the BH (Ostri ker, Weaver, Yahil, & 
|McCray||1976||Cowie, Ostriker, & Stark||1978TlB isnovatyi- 
Kogan X~Blinnikov||1980t |Krolik & London||i983| |Vitello 
1984; Wandel, Yahil, & Milgrom 1984; Ostriker, Choi, Ciotti, 



Novak, & Proga 2010 1, and because radiatio n pressure ac- 
celerates the gas and dust away from the BH (Shapiro 1973 



Ostriker, Weaver, Y ahil & McCray||1976} [Begelman( ]l985 
|Ricotti, Ostriker, & Mack||2008[ ). Indeed, the maximum ac 
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cretion rate that can be achieved in most cases is the Ed- 
dington rate. In this limit the outward acceleration of the 
gas due to Compton scattering of radiation with free elec- 
trons equals the inward gravitational acceleration. For su- 
permassive BHs (SMBHs) the picture is far more complex 
because the hot ionized region and the Bondi radius extend 
to galactic scales. Thus, it is necessary to model the evolu- 
tion of the host galaxy, with its stellar populations and the 
interstellar medium in addition to the accreting SMBH. De- 
spite the complexity of this problem, a significant amount 
of work exists in the literature on the self-regulation mech- 
anisms f or the growth of SMBHs at the centers of elliptical 
galaxies ([C iotti & Ostrik er||2001[ |Sazonov, Ostriker, Ciotti, 
& Sunyaey||2005[ |Ciotfi & Ostriker||2007| |Ciotti, Ostrikerj 
& Proga| |2009[ |Lusso & Ciotti||2011[ |Novak, Ostriker, &| 



Ciotti||201 1| |2012| l, for radiation-driven axis-symmet ric out- 
flows in active galactic nuclei (AGNs; Proga 2007; Proga 



Ostriker, & Kurosaw a 2008 1 Kurosawa, Proga, & Nagamine 
2009| |Kurosawa & Proga 2009a b ) an more generally on the 
co-evolution of galaxies and their SMBHs. The study of the 
cosmological origin and growth of SMBHs is a multi-scale 
problem often approached using large scale simulations. Typ- 
ically the limited resolution of these cosmological simulations 
precludes resolving the SMBH Bondi radius, hence, the ac- 
cretion rate is modeled using sub-grid recip es largely based 
on the Eddington-limited Bondi formula (fVolonteri & Rees 

)7U 



2005] [ Pelupessy, Di Matt eo, & Ciardi||2007[ |Greif, Jo hnson, 
Klessen, & Bromm|20081|Alvarez, Wise, & Abel 200 9t|Kim, 
Wise, Alvarez, & Abel|201 l[|Blecha, Cox, Loeb, & Hernquist| 
201 1 Blecha, Loeb, & Narayan 2012 ), sometimes with the in- 



troduction of a normalization para meter for the accretion rate 
to match observations dD i Matteo et al. 2005 ; Springel et al. 
[20051 |Di Matteo et al.|2M8j i~ 

Recently, there has been renewed interest on studies of ac- 
cretion onto in t ermediate mass BHs (I MBHs) (Milosavljevic 
eTaLl|2009a|bl |Park & Ricottil|20TT] |O1 |201 1| |Park & Ri- 



cotti 2012 1, motivated by observations of ultral uminous X-ray 
sources (DLXs) in nearby galax ies (for a review |Miller & Col 



bert|2004 van der Marel|2004| , and in the early universe, be- 



cause some simulations p redict formation of IMBHs a s Pop- 
ulation III star remnants ([Abel, Anninos, Norman, & Zhang 
TS^jBromm, Coppi & Larson||1999| |Abel Bryan, & Nor- 
man 2000; Madau & Rees 2001 ; Schneider, Ferrara, Natara- 



jan & Om ukai 2002 ) |Oh & H aiman 2002; Whal en & Fryer 
201 1 Jeon, Pawlik, Greif, Glover, Bromm, Milosavljevic, & 



Kless enU2"0"ll[ |Johnson, Whalen Fryer, & Li||2012[ |Stacy. 



Greif, & Bromm 2012) or from direct collapse of quasi-stars 
(ICarr, Bond, & Arnett||1984[ |Haehnelt, Natarajan, & Rees 
1998, Fryer, Woosley, & Heg er||2001| jB egelman , Volonteri, 
& Rees|200gnV olonteri, Lodato, & Natarajan 2008; Omukai, 
Schneider & Haiman 2008 ; Regan & Haehn eTtP009TlM ayer, 
Kazantzidis, Escala, & Callegari 2010; Johnson, Khochfar, 
Greif, & Durier|20 11 b . 



Our previous work in this series (Park & Ricotti 20 1 1 1 2012 



hereafter Paper I and II respectively), has been partially mo- 
tivated by the need to better describe the sub-grid recipe for 
accretion onto IMBH (and SMBH, with some caveats) in large 
scale cosmological simulations. We have taken the most basic 
approach to understanding the growth of IMBHs and SMBHs, 
focusing on the fundamental problem of Bondi -type accretion 
(accretion from a homogeneous and isotropic gas) modified 
by a simple feedback loop that couples the accretion rate to 
the luminosity output: L = rjMc 2 . The photons-gas cou- 



pling is due to isotropic emission of radiation from a central 
source, affecting the gas inflow through thermal and radia- 
tion pressures. We have started with idealized initial con- 
ditions, adding one physical process at a time, in order to 
understand and model analytically the modified Bondi prob- 
lem for a large parameter space of initial conditions. ID 
and 2D radiation-hydrodynamic simulations are used to ex- 
plore the large parameter space of BH masses, ambient gas 
density/temperature, radiative efficiency and spectrum of the 
emitted radiation. We have then derived a physically moti- 
vated model for the BH growth and luminosity and provide 
analytic formulae that fully characterize the model. Despite 
the simplicity of the initial conditions, we found a rich phe- 
nomenology that can be described with a simple scaling ar- 
guments. The main results are: the mean accretion rate is 
proportional to the thermal pressure of the ambient gas and is 
always less than 1% of the Bondi rate. The accretion rate is 
periodic, with period of the luminosity bursts proportional to 
the average size of the ionized hot bubble and peak accretion 
rates about 10 times the mean. Most interestingly, we have 
discovered that there are two distinct modes of oscillations 
with very different duty cycles (6% and 50%), governed by 
two different depletion processes of the gas inside the ionized 
bubble (see Paper I and II). 

This paper is the continuation of our previous work on char- 
acterizing how the Bondi-Lyttleton problem is modified by 
radiation feedback, but here we focus on the growth rate and 
luminosity of BHs in motion with respect to their surround- 
ing medium: i.e., Lyttleton accretion modified by a radiative 
feedback loop. Considering the motion of the BH is important 



for the study of stellar BH (e.g. [Wheeler & Johnson 20 IT) and 
IMBHs accreting from the ISM, both in the early univer se and 
in the local galaxies as candidates for ULXs ( Krolik, McKee, 
& Tarterj 19811 |Krolik & Kallmanj 19841 |Krolik|20041 [Ricotti 
& Ostriker||2004[ |Ricotti, Ostriker, & Gnedin|[2005l [Ricotti 
2007||Strohmayer & M ushotz ky|20091 l, and SMBHs in merg- 
ing galaxies. Both radiation feedback and BH motion are ex- 
pected to reduce the accretion rate with respect to the Bondi 
rate, but it is unknown what is the combined effect of these 
two processes. To the best of our knowledge the present study 
is the first to consider radiation feedback effects on moving 
BHs. 

This paper is organized as follows. In Section [2j we intro- 
duce basic definitions used throughout this series of papers 
and describe the numerical simulations. In Section [3] and [4j 
we show our simulation results and the analytical model de- 
scribing the set of simulations respectively. Finally, we sum- 
marize and discuss the implications of our work in Section[5] 

2. BASIC DEFINITIONS AND 
RADIATION-HYDRODYNAMIC SIMULATIONS 

For sake of consistency with the previous papers of this 
series, we define the dimensionless accretion rate A ra d = 
M /Mb, where Mb is the Bondi accretion rate for isother- 
mal gas (7 = 1). In Paper II we found that the mean accretion 
rate for a stationary BHs regulated by radiation feedback is 

(M) = min (1%T^M B , JT^Edd), (2) 

where we have defined = 7^/(10™ K). Equation Q 

is valid for density nn.oo ^ 10 5 M 2 ~ 1 cm~ 3 where a sim- 
ilar definition M n = Mbh/(10 n M©) is made. Instead, 
if ^h,oo ;$ lO 5 ^/^ 1 cm" 3 , the dimensionless accretion 
rate in the sub-Eddington regime shows a dependency with 
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FIG. 1. — Accretion rate as a function of time for a BH with Mbh = 
100 M , n H ,oo = 10 3 cm- 3 , and Too = 10 4 K. Different lines show 
simulations for BH moving at Mach number 1,2,4,8 and 10 (see labels un- 
der each line). Early phase of the simulations shows oscillatory behavior 
which becomes steady later being dominated by the motion of the gas as the 
simulations evolve. Simulations with bigger Mach number show smaller ef- 
fect of early phase oscillation which is extremely weak in simulation with 
M = 8, 10. Average accretion rate as a function of Mach number increases 
for Mach number from 1 to 4 and decreases for higher Mach numbers. Quasi- 
periodic burst of accretion is seen for the simulation with M = 4. 



the density (A ra d) = 1%^ , 4 n 5 ' , where again we define 
n m = «H,oo/10 m cm -3 . The Eddington luminosity and 
accretion rates are defined as Z^dd = 47rGMbh m p ca^, = 

3.3 x 1O 6 L M 2 , and MEdd = V 1 LEdd c ~ 2 > respectively. 

The numerical method used in this paper is the same as in 
Paper I and II. We run a set of 2D radiation-hydrodynamic 
simulations using a modified parallel versi on of the non- 
relativistic hydrodynamics code ZEUS -MP (Stone & Nor- 



|man|1992 Hayes et al. 2006 ) with our photon-co nserving UV 
and X-ray ID radiative transfer equation solver (Ricotti et al. 



|2001| |Whalen & Norman|2006] l. We include photo-heating 
photo-ionizafion, and chemistry of H/He in multi-frequency 
are manifested in the code to see how high energy UV and 
X-ray photons regulate gas accretion onto moving BHs. Ra- 
diation pressure both on electrons and H I is calculated to sim- 
ulate the effect of momentum transfer from ionizing photons 
to gas. See Paper I and Paper II for detailed description. 

Axis-symmetric geometry with respect to the azimuthal an- 
gle ((f)) is applied to all simulations. We use logarithmically 
spaced grid in the radial direction (r) and evenly spaced grid 
in the polar angle direction (0 < 9 < 7r), with BH cen- 
tered at the origin r = 0. Axis-symmetric configuration is 
necessary to simulate BHs in motion relative to ambient gas 
which is assumed to be moving parallel to the polar axes. 
In the radial direction we apply flow-in boundary conditions 
at the outer boundary for the first half of the polar angle 
(0 < 9 < 0.57r) and flow-out boundary conditions for the sec- 
ond half (0.57T < 9 < tt). At the inner boundary, in the radial 
direction, we apply flow-out boundary condition for the entire 
polar angle. Reflective boundary conditions are applied along 
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FIG. 2. — Accretion rate as a function of Mach number. Classical Bondi- 
Hoyle accretion predicts monotonic decrease of accretion rate as (A rat j) <x 
(1 + Ai 2 )~ 3 / 2 (dashed line). However, our simulations show that (A ra< j) 
decreases with M for low Mach number M < 1 with minimum (A ra< j) at 
A4 ~ 1, increases for 1 < M < Mr, and decreases again for > Mr. 
Solid line shows our model using isothermal (7 = 1) shock and D-type I- 
front jump condition for the dense shell in the upstream direction. Density 

1/2 

dependence of the mean accretion rate ((A ra d) <x n H from Paper I) is 
applied for comparison between simulations with various densities (nu,oo = 
10 2 - 10 6 cm 3 ). Dot-dashed line indicates the critical Mach number Mr 
where the (A rac i) peaks. 



the polar axis 9 = and 9 = it to satisfy the axis-symmetric 
configuration. 

We assume uniform density and constant velocity for the 
initial conditions. For the supersonic cases (VW > 1), we 
start the simulations with an assumption of fixed accretion rate 
((Arad) = 0.001) to reduce the effect of oscillation observed 
at early phases of the simulations. For our reference simu- 
lations we select a fiducial value of the radiativ e efficiency 
77 = 0.1, typical for thin disk models (Sha kura & Sunyaev] 
1973| l, BH mass Mbh = 100 M Q , and the temperature of 
the ambient gas = 10 4 K. We explore a range of Mach 
numbers up to A4 = 10, and gas densities nn.oc = 10 2 - 
10 6 cm- 3 . 



3. NUMERICAL RESULTS 
3.1. Accretion Rate as a Function of Mach Number 

Classical Bondi-Hoyle-Lyttleton accretion predicts mono- 
tonic decrease of accretion rate with increasing velocity of 
the BH with respect to the ambient gas as in Equation (fTI), 
however our simulations of moving BHs with radiative feed- 
back show a very different dependence of accretion rate as a 
function of the Mach number. 

Figure [T] shows the time evolution of accretion rate for dif- 
ferent Mach numbers .M = 1, 2, 4, 8, and 10 for gas density 
hh.oo = 10 3 cm -3 . Most simulations, except the one with 
A4 = 4, show steady accretion rates after an initial transient 
phase that is the result of out of equilibrium initial conditions. 
The time duration of the transient phase is proportional to the 
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FIG. 3. — Density and ionization fraction for simulations of a BH of mass Af^h = 100 Mg, gas density nH,, 



10 5 cm -3 , and temperature = 10 4 K 



moving at A4 = 1,2, and 4 (from left to right). Each panel shows large scale view of the H II region in cometary shape and the dense shell in the upstream 
direction for each simulation with different Mach numbers. The size of the Stromgren sphere in the downstream increases roughly linearly with increasing Mach 
number. However, the size of the Stromgren in the upstream direction remains constant as a function of Mach number. With increasing Mach number, the 
density of the shell in the upstream direction increases ((^H,sh/ n H oo) x M' 2 ) and the density behind the dense shell also increases njj i n /njj oo °c -M 2 for 
1 < M < Mr. 



crossing time scale r cr , and thus inversely proportional to the 
Mach number. The mean dimensionless accretion rate (A ra d) 
is measured when the accretion reaches a steady state while 
for the non-steady cases, such as M = 4, we calculate the 
time-averaged accretion rate. As found in Paper I, at low 
gas densities, i.e. nn,oo 5, 10 5 cm' 3 for simulations with 



sity ((A r ad) oc nVi,)- 



Mbh = 100 Mq, (A ra( j) is proportional to square-root of den- 
'h oc'- Figure 2 shows that after correcting 
(A ra d) for the aforementioned density dependence also found 
for stationary BH, the same functional form for the accretion 
rate fits all simulations with different gas densities (large sym- 
bols for rj = 0.1 and gas densities nn.oo = 10 2 -10 6 cm~ 3 ) 
and radiative efficiencies (small pentagons for a simulation 



with r\ — 0.03 and % c 



10 5 



5 ). 



The quasi-periodic oscillations of the accretion rate found 
in Paper I and Paper II for stationary BHs is still observed 
in simulations with Mach numbers M < 0.5, which main- 
tain the main characteristics of spherically symmetric accre- 
tion. This implies that introducing small systematic subsonic 
velocity to spherically symmetric accretion does not signif- 
icantly alter the qualitative description of the accretion dis- 
cussed in Paper I and Paper II. However, the average accre- 
tion rate (A ra( j) decreases steeply as a function of Mach num- 
ber, and (Arad) at M ~ 1 is roughly one order of magnitude 
smaller than for non-moving BHs, when all the other parame- 
ters are held constant. This decrease of the accretion rate with 
increasing velocity is qualitatively consistent with the predic- 
tion of Bondi-Hoyle-Lyttleton accretion. 

But the spherically symmetric accretion model fails for su- 
personic BH motion (M > 1). The shape of the H II region 
makes a transition to a well-defined axis-symmetric geome- 



try, elongated along the direction of the gas flow in the down- 
stream direction, while a bow-shaped dense shell develops in 
front of the H II region in the upstream direction, significantly 
affecting the velocity field of the gas inflow. In most simu- 
lations steady state accretion is achieved for supersonic BH 
motion, since gas is continuously supplied to the BH without 
interruption. 

Interestingly, as the BH motion becomes supersonic and a 
bow-shock and dense shell form, (A ra d) increases as a func- 
tion of Mach number. This is clearly at odds with the results 
expected from the classical Bondi-Hoyle-Lyttleton model. A 
Mach number of M ~ 1 is roughly the turning point where 
(Arad) has a minimum, while there is a maximum value of 
the Mach number M = Mr (Mr ~ 4 for = 10 4 K), 
at which (A ra d) reaches a maximum value, before staring to 
decrease with increasing BH velocity. An instability of the 
dense shell that leads to bursts of accretion rate is observed 
in some simulations in this Mach number range. This result 
will be discussed in Section |4.3| At higher Mach numbers 
(M > M r), a steady state solution is achieved once again 
since the dense shell does not form due to the high velocity of 
the gas inflow (as shown in Section 4. 1 the I-front transitions 
from D-type to R-type). For M > Mr, (A ra d) decreases 
monotonically as a function of Mach number and converges 
to the Bondi-Hoyle-Lyttleton solution shown as a dashed line 
in Figure [2] 

3.2. Structure of the Gas Flow and H II Region 

For Mach numbers 1 < M < Mr, a dense bow shock 
forms in front of the H II region in the upstream direction, 
followed by a D-type I-front (see Figure [3]). Most of the 
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FIG. 4. — Density, velocity, temperature, and H I abundance (from top to bottom) profiles along X-axis for simulations in FigurefSlfor M = 1 (solid lines), 2 
(dotted), and 4 (dashed). Left panels (note that small scale is located on the right side) show profiles of downstream while right panels show upstream profiles. 



gas inflow propagates through the bow shock without chang- 
ing direction, while a small fraction of the gas inflow in the 
outer parts of the bow-shock is re-directed farther from the 
axis of symmetry. The formation of a bow shock in the up- 
stream direction changes the gas density and velocity behind 
the shock, while the gas temperature remains relatively un- 
affected (isothermal shock). Note that the H II region has a 
cometary shape, with overall length increasing linearly with 
increasing Mach number. The size of the H II region in 
the upstream direction is not sensitive to the Mach number, 
while in the downstream direction the length of the ionized 
tail shows a linear relationship with the Mach number as 
shown in Figure [3] The upper panels in Figure [3] show the 
changes in the density structure and the H II region shape for 
M = 1,2, and 4 respectively. The lower panels show the 
vector fields over the gas density for each simulation. In the 
bottom panels, we use a logarithmic scale for the radial direc- 
tion to better show the motion of gas in the vicinity of the BH. 
The location Bondi radius calculated for the gas temperature 
inside the H II region is close to, but within the location of 
the I-front. We will see in the next section that this will allow 
us to derive simple analytical formulae for the accretion rate 
using the Bondi formula for the gas inside the H II region and 
a model for the I-front. 

For M. = 1 the size of the H II region in the downstream 
direction is roughly ~ 4 times its length in the upstream direc- 
tion. The density structure in the downstream direction is very 
complex as shown in Figure [3] The re-directed gas streams 



form high density regions and shocks. However, since most 
of the gas downstream of the BH is not accreted onto the BH, 
we will focus on understanding the upstream structures. The 
size of the H II region in the upstream direction will be dis- 



cussed in greater detail in Section 4.4 



4. ANALYTIC MODELLING 

In this section we show that properties of the gas flow 
around the moving BH can be understood using a simple ID 
model for a standing I-front in the frame of reference moving 
with the BH. There are two regimes for the I-fronts: a /?-type 
front and D-type front, determined by the BH velocity and 
density of the ambient medium. As in the previous papers we 
find that the accretion rate onto the BH can be estimated us- 
ing the Bondi-Lyttleton formula for the gas inside the I-front, 
since the effective Bondi radius lays within the ionized region. 

4. 1 . Transition from R-type to D-type I-front 

In most simulations we found that the gas flow reaches a 
steady state solution. In these cases the position of the I-front 
with respect to the BH is stationary, hence the I-front propa- 
gates with respect to the gas upstream with constant velocity 
Ubh- Although the geometry of the H II region and the bow- 
shock are clearly not planar, we can approximate the flow as 
ID (parallel to the direction of motion of the BH) near the 
I-front location in the region around the axis of symmetry of 
the problem. But we will show that accounting for deviations 
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from the planar symmetry is necessary to understand the sim- 
ulations results. 

Let's start by writing down the equations for the propa- 
gation of ID I-front through an homogeneous medium with 
density p^ . The density ratio across the I-front can be esti- 
mated by solving the mass and momentum conservation con- 
ditions Pi n v in = PooWbh = J/Ji> where J is radiation flux and 
\i = 1.27m g, assuming t he Helium becomes singly ionized 
in the front (|Spitzer|198T]): 



A, 



_ p in _ (1 + M 2 ) ± y/(l+M 2 ) 2 -4M 2 A T 



Po 



2Ai 



(3) 



Here, we have defined At = Tin/Too > 1- Due to the condi- 
tion for the density ratio in Equation Q to have real positive 
values, the Mach number must be M. < A4d, where D stands 
for dense gas, or M. > M.r where R refers to rarefied gas. 
D- and /^-critical Mach numbers are respectively: 



A T (1 - y/l - l/Ar) 
A t + 



1/A r 

Also, MdMr = 1, thus Md = M 
velocity the critical velocities are: 





I'D 



VR 



= c s , in (1 - y/l - 1/Ar) 



= c sM (1 + v/1 - 1/A T ) 



At»i 2c s 
> 2c s . 



A T >>i 



(6) 
(7) 



Only for the special case in which At — 1, we have A4d — 
A4r ~ 1 and a solution exists for any BH velocity. However, 
for all physically motivated cases with At > 1 a solution 
is not possible for BH velocities vrj < Ubh < vr. In our 
simulations if Vbh > vr ~ 2c s .i n the analytical solution that 
reproduces the data is the one with the negative sign in Equa- 
tion This solution describes a weak /?-type I-front, that 
has p; n ~ poo and supersonic motions of the gas both ahead 



and behind the I-front (in the frame of reference comoving 
with the BH). 

If Vd < c Si00 /2 < «bh < V R ~ 2c s .i n a solution is not 
possible. A shock front must precede the I-front, increasing 
the density and reducing the gas velocity below vd (£>-type 
solution). This is indeed observed in the simulations. For a 
D-type front the density past the front is always lower than the 
density upstream. Figure [4] shows that the post-shock density 
and velocity are a function of the Mach number, while the 
temperature of the shell is always T sn w T^, i.e., the shock is 
isothermal (7 = 1). The ratio between the densities at infinity 
and behind the isothermal shock is, 



Psh 

Poo 



M 2 



(8) 



that is in good agreement with the simulation results (see Fig- 
ure [5] left panel). However, for the calculation of the gas ve- 
locity in the shell our simplifying assumption of planar geom- 
etry fails. Assuming planar geometry, mass flux conservation 
in the direction of the BH motion implies p s hV s h = Poo^bh- 
Thus the Mach number in the shell Ai s h = Wsh/c s , s h is 
A4 s h = Mr 1 < 1. However, this assumption fails to repro- 
duce the simulation results and is inconsistent with obtaining 
a D-type solution for the I-front. Hence, the model must as- 
sume (as verified in the simulations) that the gas velocity in- 
side the shell has a non-zero tangential component. It follows 
that the component of the velocity along the direction of mo- 
tion of the BH is reduced due to off-axis motions that divert 
some of the gas in the direction parallel to the bow-shock. As 
illustrated in left panel in Figure : 
A4 s h ~ const, as a function of tr 



5] the simulations agree with 
le Mach number: 



M sh - AM D , with A < 1, 



(9) 



implying that v sh p sh /(p oc v hh ) = A(M/Mr) < 1. Note 
that in order to be consistent with our model, the data 
points for A4 s h are calculated from the simulation data by 
enforcing mass conservation across the I-front: A4 s h — 
PinVin/(PshC s ,sh)- This ensures that the jump conditions 
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across the I-front ignore the component of the velocity par- 
allel to the bow-shock. 

As discussed above, across the I-front, the density ratio be- 
tween the gas in the shell and in the H II region can be esti- 
mated by solving the mass and momentum conservation con- 
ditions: 



Pin Jl+M^)± y/(l+Ml 1 ) 
Psh 



4M 2 h A T 



2An 



>2M 2 D for ^4 



1. 



(10) 



For Equation ( 10 1 to have real positive values, we must have 



M s h < Md = M^ 1 (D-type solution). Clearly, if A < 1 
(i.e. M s h < Md) a D-type solution is possible for any Mach 
number M < Mr. Note that a D-type solution would not 
be possible if we had assumed ID flow geometry for which 
M s h = Mr 1 ■ ID-flow assumption would give a D-type so- 
lution only for M > Mr, that is the regime in which we 
expect an /?-type front, and is thus ruled out. In general, for 
A 2 < 1 the D-type solution consistent with the simulations 
results is the one with the negative sign in Equation dTOj, that 
has the larger relative decrease in density across the front: 
pin/psh ~ A 2 M 2 D . This solution describes a strong D-type 
front. By combining Equations ([8]), (|9| and ( 10 1, we get 



An _ 
Poo 

M 2 
2A T 



Pin Psh 
Psh Poo 

n ( M 



\Mr 



for A ~ 1, 



(11) 



or A p = A 2 {Ml Mr) 2 for A 2 < 1. The velocity ratio be- 
tween the gas inside the H II region and the BH velocity is 
Vin/v h h = (psh/pin){M s h/M) « M R /2M for A ~ 1 or 
"in/^bh = Mr/ AM for A 2 < 1. In this regime (D-type I- 
front) the model predicts that the Mach number of the gas in- 
side the H II region is a constant close to unity, as is indeed ob- 

1/2 

served in all the simulations: M- ln « (v[ n /vbh)(M/ Aji ) ~ 
1 for A ~ 1 (or M in = 2/ A for A 2 < 1). 

The symbols in the upper left panel in Figure B] show the 
ratio between the gas density at infinity and inside the dense 
shell as a function of M for a set of simulations with dif- 
ferent parameters (i.e., gas density and radiative efficiency 
as shown in the figure's legend). Our model is shown as a 
solid line. The lower left panel shows the Mach number in 
the shell M sn in units of Md, as a function of M for the 
same simulations. The right panels in FigureBJshow the den- 
sity and velocity ratios between the gas at infinity and inside 
the H II region as a function of Mach number. We measure 
the gas density and velocity within the H II region where 
the density profiles has a minimum behind the I-front. As 
discussed above, for a plane parallel I-front and shock the 
mass flux of the gas is conserved: p m v\ n / {pooVbh) = 1- 
However, due to the formation of a bow-shock when the I- 
front becomes D-type, the gas flow has a velocity compo- 
nent perpendicular to the direction of the BH motion and 
Pi-nVin/ (PooVhh) ~ AM /Mr < 1, in agreement with the 
simulation results shown in the right bottom panel of Fig- 
ure [|] In our fiducial simulations with = 10 4 K, we have 
Mr ~ 4.7 (A T ~ 6 and T in = 6 x 10 4 K). The simulation 
results confirm that a D-type I-front forms for M < Mr, 
while a transition to /?-type occurs at M ~ Mr. 

4.2. Accretion Rate 



Since the Bondi radius calculated using the density and 
velocity inside the H II region is smaller than the radius of 
the I-front, the mean accretion rate can be estimated using 
the Bondi formula for gas inside the H II region: M oc 
^hhP'^ c shX^ + -^D ~ 3 ^ 2 - A s m previous papers of this 
series we define the dimensionless accretion rate in units of 
the Bondi rate Mb oc M^poaC' 



„-3 . 

"3,00 " 



(Arad) = 



M 



Pin 

Poo 

A, 



1 



(1+^)3/2 



a 3 t /2 (i+Mir/i- 



(12) 
(13) 



We have identified two regimes corresponding to the D-type 
and R-type I-front solutions. If 1 < M < Mr, in the pre- 
vious section we found A p s» 2{M/ Mr) 2 = 2{v\ }n /vR) 2 , 
and M m ~ 1. This means that inside the H II region the 
density increases as the BH moves faster but the Mach num- 
ber remains constant (approximately transonic). This model 
explains in simple terms the increase in the accretion rate 
with increasing BH velocity observed in the simulations in 
this regime. If M > Mr, the shock front does not form 
and A p -> 1 (if M >• Mr), while the Mach number is 

1 /2 

.Mm ~ M/Ap/AJ . Hence, in dimensionless units: 



(A ra d) 



(2A T )- 5 / 2 M 2 



2 5/2 

Ml 



A^AIAt + M 2 ) 



-3/2 



M~»Mr 



I M \2 

if M < Mr 

[(Mr/2) 2 +M 2 }- 3 / 2 , 

if M > Mr. 

(14) 

The solid line in Figure|2]shows our model for the accretion 
rate onto moving BHs. The model generally is a good fit to 
the simulation results, except the model slightly overestimates 
the accretion rate around the critical Mach number Mr. 

We run a complementary set of simulations to study more 
precisely the changes of the physical properties as a function 
of Mach number, since the simulations with constant veloc- 
ities have a coarse sampling in velocity space and show an 
intrinsic scatter which is probably the result of out of equi- 
librium initial conditions. We start the simulation assuming 
sonic motion of the BH (M ~ 1) and increase gradually 
the velocity of the gas inflow at the boundary. This type of 
"wind tunnel" numerical experiments is useful to focus on the 
changes of physical properties as a function of velocity, while 
holding the other parameters fixed. The critical Mach num- 
ber Mr and the peak luminosity depends on the temperature 
ratio At as in Equation (5}. Figure 10 shows the accretion 
rate as a function of Macn number for different gas temper- 
atures at infinity = 7 x 10 3 , 10 4 , and 1.3 x 10 4 K for 
Mbh = 100 M , n H ,co = 10 5 cm -3 , and r\ = 0.1. Our 
model is in good agreement with the simulation results as 
shown in Figure 10 The peak accretion rates and the criti- 



cal Mach numbers in these cases are very close to the model 
predictions. However, the caveat is that the dense shell which 
initially forms at the beginning of the simulations does not 
change its location as the velocity of gas increases as observed 
in the simulations in which v^h was held constant. 

If we express the accretion rate as a function of the BH 
velocity in physical units, we find that the accretion rate is 
independent of the temperature of the ambient medium and 
peaks at about twice the sound speed inside the H II region: 
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FIG. 6. — Density and ionization fraction for a simulation of a BH of mass Mbh = 100 Mq, gas density 10 4 cm -3 , and temperature Too = 10 4 K moving 
at Ai = 2.7. The shell forms and gets destroyed in a quasi-oscillatory behavior. When the dense shell breaks and falls onto the BH, accretion rate shows peak 
luminosities which are much higher than the average depending on the density of the shell, and thus the Mach number Ai. 
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FIG. 7. — Simulation of a BH with mass Mbh = 100 Mq, gas density 
10 3 cm -3 , and temperature r M = 10 4 K moving at M = 3.4. As the 
Ai approaches the A4r, the dense shell in the upstream direction becomes 
unstable since the recombination time scale r rcc becomes comparable to the 
crossing time scale r cr . When the dense shells break and fall onto the BH, ac- 
cretion rate shows peak luminosities which are an order of magnitude higher 
than the average (shown as dotted line). 



v ™h X = V R — 2c Si in (see Figure 111. This is contrary to 
the case of a static BH (or moving at sub-sonic speed) for 

which (A ra d) oc 2 and hence the accretion rate M oc T^ 
is proportional to the temperature of the ambient gas. The gas 
accretion rate for a moving BH is, 

M « P°°( GMbh ) 2 x (° 7 {^z) 2 ifc s ,oo <«bh<2c s , in 

4m \ [1 + (V bh /C s ,i„) 2 r 3/2 if V hh > 2c s . in 

(15) 

The velocity for peak accretion depends only on the sound 
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FIG. 8. — Size of H II region in the upstream direction (9 = 0) {R s )e=o as 
a function of Ai. Since the recombination time scale T rec is small, {R s )g = o 
is not sensitive to the velocity of the flow. Also, it is not sensitive to the 



density of the ambient medium since (R s )g—Q oc rij 
show a good match with the model. 



-1/6 



. Simulation results 



speed inside the H II region and is vr = 50 km s 1 for 
T in = 6 x 10 4 K (c Mn = 25 km s" 1 ). We find a mild de- 
pendence of Ti n on the density of the ambient medium. As 
explored in more detail in Paper I, Tj n depends on the hard- 
ness of the spectrum emitted by the BHs, the gas metallicity 
and Compton cooling/heating. 

Our model has a peak accretion rate of about 70% of the 
Bondi rate in a gas with temperature T; n and density pco, 
in good agreement with our ''wind tunnel" type simulations 
(however, as noted before, the simulations with constant BH 
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FIG. 9. — Size of H II region in the downstream direction as a function 
of M. Maximum size of H II region and (R s )g =n are shown as squares 
and triangles respectively. Our model shown as a dashed line approximately 
matches with the maximum size of the H II region of the downstream. The 
size of the H II region in the direction 8 = tt is smaller than the maximum 
size of the H II region due to the concentrated gas along the line of 6 = 7r. 

Mach number produce smaller accretion rates at peak accre- 
tion: about a factor of 5 smaller than the model prediction for 
the peak value). Interestingly, in our model and simulations a 
BH moving at 20-50 km s _1 with respect to a gas with tem- 
perature ;$ 10 4 K, has a faster growth rate and accretion 
luminosity than if it was at rest (or moving at subsonic speed). 
Figure 1 1 is an extrapolation of our results to lower tempera- 



ture regime showing accretion rates as a function of BH veloc- 
ity Vbh for gas temperatures T^ = 10 4 , 3 x 10 3 , and 10 3 K. In 
all cases, if Too < 10 4 K, we find that the peak accretion rate 
at «bh ~ 2c Si i n is larger than the corresponding accretion rate 
onto a stationary BH. Clearly this is an important result be- 
cause significant BH accretion is only possible when the BH 
is in a dense medium, for instance a molecular cloud or the 
cold neutral medium (CNM) in a galaxy, that generally have 
low temperatures. In addition, the colder the gas the smaller 
the BH velocity needs to be to achieve the supersonic speed 
that leads to significant increase of the BH accretion rate. 

In the next section we will discuss another important re- 
sult: in the regime when the I-front is D-type, instabilities of 
the dense and thin shell behind the bow shock may lead to 
its fragmentation and hence produce periodic oscillations of 
the accretion rate and luminosity of the BH. This effect can 
further increase the peak accretion luminosity of the BH by 
roughly a factor of ten (i.e. p s h/Poo ~ M 2 ) with respect to 
the mean values estimated above in equation dT5b. 



4.3. Stability of Bow Shock and Periodic Oscillations of the 
Luminosity 

As discussed in the previous section, the average accre- 
tion rate (A ra d) increases with increasing Mach number if 
1 < M < Mr. For the lower values in this Mach num- 
ber range, all simulations approach a steady state accretion 
rate. Interestingly, as the Mach number approaches Mr, 
the bow-shock in simulations with ambient gas densities in 
the range tih.oo = 10 3 -10 4 cm -3 , becomes unstable pro- 
ducing intermittent bursts of accretion due to a cyclic for- 
mation/destructi on of a dense shell in the upstream direction 
(see |Whalen & Norman] 201 1} . Figure [6] shows time evo- 
lution of a simulation with Mhh = 100 M™, gas density 
^h,oo = 10 4 cm -3 , temperature Too = 10 K and Mach 
number M = 2.7. As seen in Figure [3] the ionizing radia- 
tion creates a "cometary-shaped" H II region around the BH. 
In the early stages of the simulation, the gas flow remains 
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FIG. 10. — Accretion rates as a function of Mach number for simulations 



with M bh = 100 M Q , n H ,, 



10 cm , and r\ = 0.1. Different lines 



and symbols show the models and simulations respectively for temperatures 
at infinity Too = 7000 K (long-dashed line, open hexagons), 10000 K (solid 
line, open/solid squares), and 13000 K (short-dashed line, triangles). Solid 
squares show the simulation with higher BH mass Af^h = 1000 Mq and 



lower ambient gas density nu,< 



10 cm d which shows the identical 



results with the simulation with the same ambient gas temperature T x 
10000 K. 



relatively steady. However, with time some instabilities start 
growing leading to the fragmentation of the shell. As a result 
the ionizing radiation and the hot gas inside the H II region 
are no longer contained by the dense shell downstream of the 
bow-shock and an explosion takes place. Fragments from the 
broken dense shell fall onto the BH significantly increasing 
the accretion rate, thus creating more ionizing photons that 
blow out further the thinner parts of the shell. However, after 
a time delay from the burst roughly estimated as the H II re- 
gion sound crossing time, the dense shell re-forms, resetting 
the initial conditions for the next burst cycle. 

Since the dense shell is located outside the inner Bondi ra- 
dius, the gravitational acceleration on the shell is relatively 
small and hence the time scale for Rayleigh-Taylor instabil- 
ity is long. In addition the radiation has a stabilizing effect 
as it tends to smooth out the growth of linear perturbations 
on small scales, that have the faster grow rate. However, 
when t>bh > c Si i n (i.e., M ~ Mr/2) the ram pressure be- 
comes comparable to the thermal pressure inside the H II re- 
gion and the I-front is pushed closer to the BH. In this case 
the increased gravitational acceleration on the shell and the 
increased sharpness of the pressure gradient seem to trigger 
the growth of instabilities. Also the column density of dense 
shell is roughly constant as a function of the BH velocity, 
hence because p s h oc M 2 the thickness of the dens e shell is 
oc M ~ 2 . A thin shell is well know n to be unstable ( Vishniac 
1994; Whalen & Norman 2008b a). 



As shown in Figure|5J in simulations at intermediate densi- 
ties (riH.oo = 10 3 -10^;m~ 3 ) the Mach number in the dense 
shell is close but slightly larger than M d while for higher 
densities (10 5 cm -3 ) is always less than Md- This implies 
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that D-type solution for 10 5 cm" 3 exists for all Ai since 
A4 a h < M-d while a transition between D- and /?-type is 
expected for 7Jh.oo = 10 3 -10 4 cm -3 . As shown in Figure|7j 
for simulations with nn,oo = 10 3 cm -3 the interval between 
the bursts of accretion is fairly constant At ~ 3000 years. 

In Paper I and Paper II we found a linear relationship be- 
tween the average size of the Stromgren spheres (R s ) and 
the period between oscillations for stationary BHs, which is 
of a few 10 3 years for M hh = 100 M Q and r\ = 0.1. We 
have shown that this time scale is proportional to the sound 
crossing time of the H II region and have provided analyti- 
cal relationships for the period. The same argument can be 
applied for interpreting the period between intermittent bursts 
for moving BHs and in the next section we provide a model 
that describes the the size and shape of the H II region around 
moving BHs. 

4.4. Size o/H II Region in the Up/Downstream 

The gas inflow in the direction of the polar axis 9 — (up- 
stream) and 9 = it (downstream) can be approximated by a 
ID flow with the H II region being supplied with neutral gas 
with constant velocity Vbh- The total number of ionizing pho- 
tons emitted by the BH must equal the number of hydrogen 
recombinations inside a radius (R s ) in addition to the flow 
trough the I-front of neutral gas of density n and velocity v^: 



9 a T( 



-nl + 4Tr(R s ) 2 g nv hh cos(9), (16) 



where AT ion is the number of emitted ionizing photons, being 
directly related to the luminosity of the BHs (which is a func- 
tion of Mach number). When the magnitudes of the two terms 
on the right side of the Equation ( 16 1 are compared, at 6 = 
the first term is dominant over the second term due to the BH 
motion. The electron number density inside the H II region 
for 1 < M < M R is 



•^e^H.in 



M 2 



(2A T ) 5 / 2 



(17) 



Since (A ra d) oc n 



1/2 
H,oo 



(for n H .oo < 10 cm 3 and M bri 



100 M©), the average size of the H II region in the upstream 
direction is 

{R s )e= <x v 1/3 n-^. (18) 



Figure[8]shows the size of H II region at 9 = as a function of 
Mach number for simulations with various densities (7Jh,oo = 
10 2 — 10 6 cm" 3 ) and radiative efficiencies (r/ = 0.1,0.03). 
The model is a good fit to the simulation results. 

We model the size of the H II region in the downstream 
direction in a similar manner. In the downstream direction 



9 = 7r the second term on the right hand side of Equation ( 16 1 
dominates over the first term: 



47r<E s )^ =7r n H ,: 



(19) 



where n>H,i n can t> e calculated simply using pressure equilib 
rium condition nn in = n n 

in the downstream direction for 1 < M. < A4r is 



A T 1 . The size of the H II region 



{Rs)e=^^n 1 ^ 00 {l + M 2 ) l '\ 



(20) 



where (i? s )g =7r is approximately proportional to the Mach 
number. Figure [9] shows that the model reproduces the length 
of the tail of the H II region in the downstream direction in 




v /c 

00/ s, in 

FIG. 1 1 . — Accretion luminosity normalized by Eddington rate as a function 
of the BH velocity D^h in units of the sound speed inside the H II region 
for a BH with M bh = 100 Mq, n HiCO = 10 4 cm" 3 , and T) = 0.1. 
Different line types are used for models with different temperatures Too = 
10 4 (solid), 3 X 10 3 (short dashed), and 10 3 K (long dashed). Note that 
the accretion rates peak where ~ 2c s .i n regardless of Too. Dot-dashed 
lines for each temperature Too show possible enhanced luminosities due to 
the accretion of high density gas from the broken dense shell remnants (if 
p sh oc M 2 ). 



simulations for nn,oo = 10 5 cm" 3 . The length of the H II re- 
gion at 9 = 7r is shorter than the model prediction because of 
the gravitational focusing of the gas by the BH that increases 
the gas density on the axis of symmetry behind the BH. 

5. SUMMARY AND CONCLUSIONS 

In this third paper of a series on radiation-regulated accre- 
tion onto BHs, we have focused on the effect of the BH mo- 
tion relative to the surrounding gas, i.e., the Hoyle-Lyttleton 
problem modified by the effects that photo-heating and radi- 
ation pressure from the radiation emitted near the BH have 
on the accretion flow. As in previous papers of this series 
we have used radiation-hydrodynamic simulations to explore 
a large parameter space of initial conditions to inform us on 
how to formulate an analytical model that reproduces the sim- 
ulation results. The followings are our key findings. 

• The quasi-periodic oscillation of the accretion rate ob- 
served in simulations of non-moving BHs are only ob- 
served for subsonic motions of the BH, while the ac- 
cretion rate becomes steady at supersonic velocities (in 
most cases). 

• In the supersonic regime we observe an axis-symmetric 
gas flow and a "cometary shaped" H II region with tail 
length proportional to the BH velocity «bh- F° r ve " 
locities c s oo < u^h < V R ~ 50 km s" 1 the ionization 
front becomes D-type: a bow-shock and a dense shell 
develop in front of the H II region in the upstream di- 
rection. For Vbh > vr the bow-shock disappears and 
the I-front becomes R-type. 
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• For subsonic motion of the BH we find that the accre- 
tion rate onto the BH decreases with increasing BH ve- 
locity. However, contrary to naive expectations, the ac- 
cretion rate increases with increasing BH velocity in 
the regime c soo < < 50 km s _1 , when the I- 
front is D-type. The accretion rate peaks at BH veloc- 
ity i>bh ~ 50 kms" 1 before it starts decreasing with 
increasing BH velocity, converging to the well known 
Hoyle-Lyttleton solution without radiation feedback. 

• Based on the simulation results we formulate a simple 
analytical model of the problem based on modeling the 
jump conditions across the I-front. The transition of the 
I-front from R-type to D-type happens at Vhh = 2c Si ; n , 
where c s j n is the sound speed inside the H II region. An 
isothermal jump conditions for the bow-shock repro- 
duce fairly well the gas flow properties in simulations 
with ambient gas density njj.oo ^ 10 2 cm~ 3 . Because 
the inner Bondi radius (for the gas inside the H II re- 
gion) is comparable but smaller than the I-front radius, 
the accretion rate is well reproduced in our analytical 
model assuming Bondi-type accretion from the ionized 
gas downstream of the bow-shock. The BH moves sub- 
sonically or transonically with respect to the gas down- 
stream of the bow-shock. In this regime the accretion 
rate increases with BH velocity because the density of 
the ionized gas inside the H II region increases with in- 
creasing velocity reaching n# ~ iih,oo at i>bh = 2c Si ; n - 

• Simulations of BH accreting from a high-density 
medium (nn.oo = 10 5 — 10 6 cm -3 ) show steady ac- 
cretion rate for all Mach numbers. However, at inter- 
mediate densities (nn.oo = 10 3 — 10 4 cm~ 3 ) we find 
intermittent bursts of accretion rate in the Mach num- 
ber range 2.5 < M. < M.r. The oscillatory behavior 
of accretion rate is due to the development of instabili- 



ties in the thin shell in the upstream direction that cause 
its cyclic fragmentation and re-formation. In lower den- 
sity regime, nn.oo ^ 10 2 cm" 3 , the post-shock density 
is lower than expected assuming isothermal shock jump 
conditions; as a result the dense shell is thicker and less 
prone to instabilities. For t>bh > Vr the dense shell 
never forms, providing steady accretion rates. 

• Contrary to the case of non-moving BHs in which 
M oc Too, the accretion rate onto supersonic BHs is in- 
dependent of the temperature of the ambient medium. 
It follows that if Too < 10 4 K the accretion rate onto a 
BH moving with velocity Vhh ~ 50 krns" 1 is about 
5(10 4 K/Tqo) times larger than the accretion rate at 
Vhh = 0. Hence, the growth rate and the mean accre- 
tion luminosity of BHs moving supersonically can be 
significantly larger than that of non-moving BHs with 
the same mass and accreting from the same medium. 

The new results reported in this paper may have significant 
repercussions in modeling and understanding the growth of 
stellar BHs in the early universe, the build up of an early 
X-ray background and provide important clues for explain- 
ing ultra-luminous X-ray sources (ULXs; Krolik, McKee, & 



|Tarter||198H IKrolik & Kallmanj|1984l |Krolik||2004| IRicotti 



& Ostriker 2004; Ricotti, Ostriker, & Gnedin 2005; Ricotti 
|2007[|StrohmIy er & Mushotz Eyffi09) . 
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